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Abstract. Nonlinear Analysis B: in press, 2004. 

The Swift-Hohenberg fluid convection system with both local and non- 
local nonlinearities under the influence of white noise is studied. The ob- 
jective is to understand the difference in the dynamical behavior in both 
local and nonlocal cases. It is proved that when sufficiently many of its 
Fourier modes are forced, the system has a unique invariant measure, or 
equivalently, the dynamics is ergodic. Moreover, it is found that the num- 
ber of modes to be stochastically excited for ensuring the ergodicity in the 
local Swift-Hohenberg system depends only on the Rayleigh number (i.e., 
it does not even depend on the random term itself), while this number for 
the nonlocal Swift-Hohenberg system relies additionally on the bound of 
the kernel in the nonlocal interaction (integral) term, and on the random 
term itself. 
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1. Introduction 

Density gradient- driven fluid convection arises in geophysical fluid flows in the 
atmosphere, oceans and the earth's mantle. The Rayleigh-Benard convection is 
a prototypical model for fluid convection, aiming at predicting spatio-temporal 
convection patterns. The mathematical model for the Rayleigh-Benard convec- 
tion involves the Navier-Stokes equations coupled with the transport equation 
for temperature. When the Rayleigh number is near the onset of the convec- 
tion, the Rayleigh-Benard convection model may be approximately reduced to 
an amplitude or order parameter equation, as derived by Swift-Hohenberg [34j. 

In the literature, most works fe.g.. |17t I14[ I28j) on the Swift-Hohenberg model 
deal with the following evolution equation for order parameter u{x,t), which 
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is localized version of the model originally derived by Swift-Hohenberg [HI] , 



where q measures the difference of the Rayleigh number from its critical onset 
value and the cubic term is an yet "approximation" of a nonlocal integral 
term in the original Swift-Hohenberg model j33]. 

Roberts ([321 EH]) re-examined the rationale for using the Swift-Hohenberg 
model as a reliable simplified model of the spatial pattern evolution in fluid 
convection. He argued that, although the localization approximation used in 
(1) makes some sense, the approximation is deficient in describing some basic 
features of such systems, and devised via symmetry argument the following 
modified Swift-Hohenberg equation with nonlocal interactions 



where G*v? is a. spatial convolution integral and G(-) is a given radially sym- 
metric positive (or nonnegative) function. We call this function G the kernel 
for the nonlocal term. In fact, nonlocal integral terms often appear naturally in 
amplitude equation models for nonequilibrium systems; see, e.g., P^IT^ ISII^. 

Our goal in this paper is to examine the above local and nonlocal Swift- 
Hohenberg models, by investigating the dynamical difference of both models 
under random impact as well as under nonlocal interactions. 

Fluid systems are often subject to random environmental influences. On the 
one hand, there is a growing recognition of a role for the inclusion of stochas- 
tic effects in the modeling of complex systems. Randomness can have delicate 
impact on the overall evolution of such systems, for example, noise-induced 
phase transition or stochastic bifurcation stochastic resonance jTHj, and 
noise-induced pattern formation |2] . Taking stochastic effects into account 
is of central importance for the development of mathematical models of such 
complex phenomena in engineering and science. Macroscopic models in the 
form of partial differential equations for these systems contain such random- 
ness as stochastic forcing, uncertain parameters, random sources or inputs, 
and random boundary conditions. Stochastic partial differential equations 
(SPDEs) are appropriate models for randomly influenced spatially extended 
systems. On the other hand, the inclusion of such effects has led to interesting 
new mathematical problems at the interface of probability and partial differ- 
ential equations. 

In this paper, we consider the Swift-Hohenberg model, taking stochastic forc- 
ing as well as nonlocal interactions, into account. 
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First, we consider the local stochastic Swift-Honenberg system (LSSH) 

(1) ut = gu - {1 + d^^fu - + F{x,t), u{0) = Uq, 

under the periodic boundary condition on with period 27r, in which the sto- 
chastic force is taken to be of the form 

N' 

(2) F{x,t) = dtWix,t) = dtY,bieiix)wi{t), N' < oo, 

i=l 

where {cj} is an orthonormal basis of the Hilbert space H = L^^^^O, 2n) formed 
by the eigenvectors of the operator A = —(1 + dxx)"^, {wi{t),t > 0} are inde- 
pendent standard Wiener processes, and the real coefficients 6j > satisfies 

bi^O , l<i<N<N' 

for some sufficiently large A^. 

Then, we study the nonlocal stochastic Swift-Hohenberg system (NLSSH) 

(3) ut = gu — {1 + Oxx^u — uG * + F{x,t), u{0) = uq, 

with a positive kernel (i.e., G > 0) and a special nonnegative kernel (i.e., 
G > 0), respectively, under the periodic boundary condition on /, in which 

G*u^= / G{x-Ou'^{^,t)d^, 
Jo 

and the stochastic force is taken to be of the form 

N' 

(4) F{x,t) = dtW{x,t) = dtY,biei{x)wi{t), N' < oo. 

i=l 

The diagnostic tool we choose to compare the dynamical differences between 
the local and nonlocal stochastic Swift-Hohenberg systems is ergodicity. Namely, 
we are interested in the ergodicity of the Markov solution process in H gen- 
erated by the stochastic Swift-Hohenberg equation which forms a random 
dynamical system (RDS). More precisely, we first prove that if sufficiently 
many of its Fourier modes are forced, the LSSH system has a unique invariant 
measure, or equivalently, the dynamics is ergodic in the phase space, and that 
all solutions converge exponentially fast in distribution to this unique measure 
(Theorem 15.21 and Corollary 15. 3p . Our results show that under certain 
conditions the RDS generated by (0) has a global exponentially attracting 
fixed point in the sense of distribution. And this fixed point is the unique 
invariant measure which is supported on the random attractor. Namely, the 
LSSH system (|T)) is ergodic. To prove our results, we decompose the LSSH 
system into a low mode part with a finite dimension and a high mode 
part with an infinite dimension. In the low mode part, we use the maximal 
coupling approach, which is used in to prove that the low mode part can 
be coupled. The high mode part can be enslaved, and then the ergodic result 
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is obtained; see also [ini EI! • 

We verify that the dynamics of the nonlocal stochastic Swift-Hohenberg equa- 
tion (jS)) is also ergodic, provided sufficiently many of its Fourier modes are 
forced (Theorems 16.41 and 16. 5p . However, our interest here is to investi- 
gate the difference in the conditions ensuring the ergodicity of the LSSH and 
NLSSH systems. We find out that (see discussions in Sec. 7): 

(i) For the local Swift-Hohenberg system, the estimated number of Fourier 
modes to be randomly excited for ensuring ergodicity depends only on the 
parameter g, which measures the difference of the Rayleigh number from its 
critical onset value. Note that this number does not depend on the random 
forcing term; 

(ii) For the nonlocal Swift-Hohenberg system with positive kernel G in the 
nonlocal nonlinearity, the estimated number of Fourier modes to be randomly 
excited for ensuring ergodicity depends additionally on the upper and lower 
bounds of the positive kernel, and on the random term itself. 

Recently, there have been a number of papers on ergodicity of stochastically- 
forced partial differential equations (SPDEs) and, in particular, 2D Navier- 
Stokes equations. The stochastic force of the SPDEs may be in the form of |21 
above or in the following form (so-called kick-force) 

oo oo 

(5) F{x, t)= J2 f'^^^t - ' = E ^^•^i'^^J- 

fc=— oo j=l 

where hj > are some constants such that ^ 6^ < oo and {ijk} are inde- 
pendent random variables with fc-independent distributions, or the white (in 
time) force of the form A coupling approach and the so-called Kantorovich 
functional (e.g., [20112211211) have been developed to show exponential conver- 
gence to the unique invariant measure for stochastic Novier-Stokes equations. 
Similar results were obtained in [201 • For the white noise forced PDE, ^JJ and 
[201 obtained the ergodicity of the stochastic forcing Navier-Stokes equation 
for the case when the random force is singular in x. For the two-dimensional 
Navier-Stokes equation with periodic boundary condition and random forc- 
ing, [9^ proved uniqueness of the stationary measure under the condition that 
all "deterministic modes" are forced, by studying the Gibbsian dynamics of 
the low modes. The idea in [01 has been used in other papers to get explicit 
results. Moreover, ^01 obtained the exponential mixing property of the sto- 
chastic Ginzburg-Landau equation by the so-called binding construction. We 
also mention that [211 developed the idea in |221 to get the rates of the 
measure converge to the unique invariant measure for 2D Navier-Stokes equa- 
tion with a white noise to all Fourier modes. 
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In the present paper, the hnear part — (1 + dxxYu in the Swift-Hohenberg 
models is not dissipative, i.e., it has positive eigenvalues. Similar situation has 
also been considered in tl6j, for example. In this case, the coupling approach 
can still be used to investigate ergodicity. In fact, we can directly study the 
maximal coupling solutions of the stochastic Swift-Hohenberg system, without 
constructing the Kantorovich functional; see fHl HZj • Then the key fact that 
the orbit in set S{m, k) comes from R{m) U S{m, m) helps us to prove the 
main result without using the exponential decay of a Kantorovich type func- 
tional, for details see Sec. 4 and Sec. 5. For simplicity of presentation, in this 
paper we only work with the periodic boundary condition and only consider 
the case when only finite modes are forced, i.e, A^' < oo. However, we empha- 
size that our argument applies in principle to the case when the other more 
physical boundary conditions are posed and when all the modes are subject 
to random excitation. Moreover, our techniques also applies to more real- 
istic two-dimensional Swift-Hohenberg model, i.e., (1 + d^x)'^ is replaced by 



We remark that a systematic treatment for random dynamical systems (in- 
cluding ergodicity) is emerging; see, e.g., PQ [7| El EHl ES] ■ 

The remainder of this paper is organized as follow. In Sec. 2, we decompose 
the local system LSSH into the low mode part and the high mode part. Sec. 3 
is devoted to some energy estimates of the solution. The coupling construction 
and a coupling result are given in Sec. 4. Then we give the ergodic result in Sec. 
5. In Sec. 6, we study the ergodicity of nonlocal system NLSSH with positive 
or nonnegetive kernels of nonlocal nonlinearity, respectively, by the techniques 
developed in the previous sections. Finally, we discuss the difference between 
LSSH system and NLSSH system by analysing the minimal conditions ensuring 
the ergodicity in Sec. 7. 



Let A = —(1 + dxxY- It is well-known that A generates a compact analytic 
semigroup e*"^ in C° j,. Since 



(1 + dxx + dyyf. 



2. Decomposition of the Local System 




N 



(9) 




i=l 



the stochastic convolution J^- 



* e^(* '^)dW{x,t) is continuous for (x, t); see jH]. 
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For convenience, we assume A^' < oo. Let (-, ■) be the scalar product of 
H = Ll^riO, 2n), I ■ I the norm, and V = H^nH. Let 

Hi = span{ei, < i < N}, = span{ei, i > N}, 

which are two subspaces of H. We will call Hi the low mode space and H^ the 
high mode space. Clearly, H = Hi^Hh- Let Pi and Ph be the orthogonal 
projections onto Hi and Hh, respectively. Denoting u{t) = {l{t),h{t)), then 
(jB))-® is equivalent to the following equations on H 

(10) / = PiAl + gl-Pi{l + hf + PiF{x, t) , /(O) = /o; 

(11) h = PhAh + Qh- Phil + hf , h{0) = ho, 

where Iq = Piu{0) and /iq = Phu{0)- For any given T > 0, since the noise is 
additive we can treat the equation pathwise. Then the deterministic theory 
implies that (jHl)-® admits a mild solution with transition probability P{t, x, F) 
for F G S{H), the Borel a-algebra in H. Here, we treat ®-(jHI) in another way 
that for Iq G Hi, ho G Hh, (fTUjl - lfTH) has a unique solution (/, h) E X = 
C{[0, T],H), where I G C(0, T; Hi) and h G C(0, T; Hh). Notice that, for given 
l{t) we can solve (10) with ho G Hh- This solution is written as $(Z, ho). Much 
more such decomposition is discussed in [TUj . 

Let li = I, hi = h for t G [0,T], 1^ and h^ be the solutions on [0, kT], lm,k and 
hm,k be the solutions on [mT, kT] for < m < fc. Similar for and Um,k. 
For (jlOp and (jllj) . we have the following proposition, which displays the main 
idea of the paper. 

Proposition 2.1. For any Uq, Mq G H, if 1]^^^ — ^mk ^^^^ there exists a 
constant A > such that 

\u^-u'^\ < e-^^'-"'^'>\u\mT) - u\mT)\ , t G [mT,kT]. 

Proof. Since ll^k = I'^k^ only consider the high mode. Let p = h^ — h^ . 
Then p satisfies 

p = Ap + Qp-Ph{{u')'~{u')'). 
Taking the scalar product of the above equation with p in Hh and notice that 

{Ph{{u')^-{u')^),p)>Q, 

we have 

dt\p\'^ < (otv + q)\p\^, 
where otv is the Nth eigenvalue of A. If we choose so large that + ^ < 0, 
then 

|p|2 < e-^(*~'"^)|p(mT)| , t G [mT,kT], 
for some constant A. This completes the proof. □ 

Remark 2.2. This proposition implies that, if we want to get the distributions 
of the LSSH system converged to the unique invariant measure, it is enough 
to prove the following: for any two solutions and , the probability of the 
event that those low modes coincide is exponentially tends to 1, as t ^ oo. 
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But, in fact, we do not need that all solutions satisfy the above proposition. By 
the coupling result in Lemma 4-1 below, we just need that the solutions that 
do not grow fast satisfy the proposition. This will also be used for the NLSSH 
system. 

3. Some Estimates 

In this section, we derive some estimates of the solutions for the LSSH sys- 
tem, which will be used in the following. We will work on a probability space 
{fl, P) generated by We associate fl with the canonical space gener- 

ated by all dwi{t). Denote as the associated a-algebra generated by 
with P the probability measure. Expectations with respect to P will be de- 
noted by E. 

The following lemma describes the growth rate of 

Lemma 3.1. There is a positive constant Ci, such that for any r > 0, 

'{l«Wr + j \u{s)\'^ds<\u{0)\'^ + Cit + r ,fort>0^ >l-e-'^. 



P 



Proof. Applying Ito's formula to one yields 

lj^Ht)\^ = {^%u) + {qu,u) - {u^,u) + ^Bo + {F{x,t),u). 
Note that 

{Au + gu, u) = (—(1 — g)u — 2dxxU — d^^u, u) 

< -(1 - g)\u\^ + \dxxu\^ + \u\^ - \dxxu\^ 

— qW\ 



|2 



and 

\u\ < (27r) 4 |m|^4. 

Then we have 



! \u{s)\''ds < |m(0)|2+ / [2(^+ ^|«|4 + So]cis + 2 / {F,u) 
Jo Jo 2 2n Jq 

= \um' + j\2{g + hl^ + \)\u\' - ^1^1' + Bo]ds + M, 
-(M),/2 - ( f 2hlM^ds - (M),/2) 



where femax = maxj bi,Mt — 2 {F, u) and {M)t is the quadratic variation. 
Note that 

{M),/2^2j2bl [ ulds<2bl^f \u\^ds. 

Jo Jo 
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That is, 



/ |u(s)pds< |M(0)|2 + Cit + Mt- (M)i/2, 





where Ci = max^{2{g + 6^^^ + ^x"^ - ^x'^ + Bq}. 
Then the classical supermartingale inequality implies 

p||M(t)P + ^ \u{s)\^ds < \u{0)\^ + Cit + r 

> p|Mt-(M)t/2<r} 

> p|exp(Mt- (M)4/2) <e^} 

This completes the proof. □ 
We also need the following estimation of the mean value of the solution. 

Lemma 3.2. For any t > 0, we have 

E|M(t)|2 < e-"*E|MoP + i?, 
where a and R are some positive constants. 

Proof. Applying Ito formula to |Mp and taking the mean value, we find 

E|Mp + a / E\u\^ds < E|moP+ / [i2g - a)E\u\^ - —E\u\^ + Bo]ds 
Jo Jo 27r 

< E\uo\'^ + Rat, 

where a is some positive number and Ra = maxx{{2g — a)x'^ — ^x^ + Bq}. 
Then Gronwall inequality yields the results. □ 

4. Coupling Approach for Ergodicity 

We start this section with some notations and terminology. Let if be a sep- 
arable Hilbert space with a-algebra B{H) and let J^{H) be the space of sigh 
measures with bounded variation. We denote by V{H) the set of probability 
measures fi G M.{H). For /i G V{H) we define the variation norm as 

\\lj\\var= sup \fl{A)\. 

AeB{H) 

If fj,i, fi2 £ '^{H) are absolutely continuous with respect to a fixed Borel 
measure m, then we have 

I l/il - /i2 1 kar = ^ / \Pl{u) - P2{u)\dm{u), 

where pj, i = 1, 2, is the density of /ij with respect to m. For more see pOj . 
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The coupling is a well-known effective tool for studying finite-dimensional 
Markov chains 12^. To our knowledge, j2H] is the first paper using a cou- 
pling approach for invariant measures of the stochastic heat equation. The 
main idea of coupling is to construct two random variables and ^2, for two 
measures /ii and /i2 that we are concerned with, and study the two measures 
through constructing the two random variables. In this section, we will use 
the maximal coupling, namely, we construct ^1 and ^2 such that 

About the maximal coupling, we refer to [20] and [21]. 

In this section, we always use the maximal coupling {ll, If) for the distributions 
of {Piu^.Piu^) on C{0,T; Hi). By the discussion in Sec. 2, we can define the 
coupling solution u\t) as 

u\t)={ll{t),^llh\{k-l)T))y te[ik-l)T,kT], k>l, 1 = 1,2 
for (5)-(7) on the product space JF'^, P^) x X'^, where 

^]'= = ^]x•■■x^],J^^'= = J^^x■■■xJ^^,p'= = Px■■■xP. 

and = X X ■ ■ ■ X X. From Theorem 4.2 in [20], it follows that is mea- 
surable coupling solutions of (0)-®- Notice that, the coupling solutions are 
not the solutions of ©-(El), but they have the same distributions with the so- 
lutions of ©-(0), since we just change F{x,t) as another noise with the same 
distribution. Then, we will call the coupling solutions the weak solutions as in 
PT] . Obviously, the results in Sec. 2 still hold true if u{t) is a weak solution 
of 

Now, for those weak solutions we define some function spaces which play a 
role in our approach. This arises from Remark 2.2. These spaces are also 
introduced in [21]. For any given real number T > and integer number 
> 1, we define So(m,k) as the set of functions {u^ (t) , (t)) , t G [mT,kT], 
< m < k, such that 

(12) u' e L'^{0,kT;V)nC{0,kT;H), 

(13) Piu\t) = Piu\t) , \u\mT)\ V \u\mT)\ < D, 

(14) S'{t,mT) <r + {Ci + l){t-mT) , mT < t < kT , i = 1,2, 

where a\/ b = max(a, 6) for a, 6 e M and £'-{t,mT) = |M'(t)P + J^j, |-u*(s)p(is. 
The set So{m, k) is called the coupled set. Then the set 

R[k) = (C(0, kT; H) n L2(0, kT; V)) \ U^^^^Soim, k) 

should be the exponentially small set by Proposition 2.1. 
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Let 

S{m, k) = So{m, k) \ S{m — l,k), < m < k 

where S{—l,k) is the null set. In the following, we always assume |mqP and 
I^Qp have the finite mean value. We observe that the orbit in S{m,k — 1) 
enters either S{m, k) or R{k) U S{k, k) for < m < — 1, and the orbit in 
S{m, k) comes from the set R{m) U S{m,m); see Fig.l. 



S(0,1)— S(0,2) 
8,(1,2) 




S(2,2) 
U R(2) 



S(m-l,m 




S(m,m 
U R(m 



S(0,m+1)- 
S(l,m+1)- 

-S(m-l,m+l) 
m,m+l) — 



S(m+l,m+l) 
- U R(m+1) 




S(k-2,k-l 



S(k-'i;k-l 
U R(k-l) 



S(k,k) 
U R(k) 



Figure 1: Orbit transition 

Then, what needs to be proved is that the event that and are in 
R{k) U S{k, k) has exponentially small probability. Clearly, this demands that 
the events that the orbits in S{m, k — 1) enter R{k) U S{k, k) or do not enter 
S{'m, k) have exponentially small probability. 



Lemma 4.1. Fix T > large enough. Let Ai and A2 he the distributions of 
ll^f, and link' respectively, for {u^^u^) G S{m,k — 1), < m < A; — 1. // (OP 
holds and N is large enough, then 

I |Al — A2\\var S 6 ' 

for some positive constant 7. 

Proof. Let and P be the maximal coupling for Ai and A2, respectively. There 
are two cases to be distinguished for t E [0, T]. 

Case l.£\{k - 1)T + t, mT) < r + (Ci + 1)((A; - 1)T + t - mT), i = l,2. 
Let 

/^ = W(mT), h}^ = Pnu\mT), i = l,2. 
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Then Iq = = I- Consider the LSSH system in low mode space 

Let B{1, hi, hi) = Piil + h^ - Piil + h^ f. Then 

mXX)? < sup \{P,(l + h'f-Pi{l + h')',y)\' 

y&L'j and \y\=l 
< |(/ + /,1)3_(/ + /,2)3|2 

Proposition 2.1 tells us that 

Then, there exists some positive constant A', such that 

since at most increases polynomially. Thus 

(15) / \B{lXX)?dt<KQe-^^^-'^-'^'^ 
Jo 

for some constant Kq > 0. Now we apply the Girsanov's formula to Aj. Let 

P{t,u;) = h~'B{l,hl,hl), 

where b is the N xN diagonal matrix with diagonal elements bj, j = 1, ■ ■ ■ ,N. 
Then, Girsanov's formula yields 

(16) Xi{dl) = e^^^^X2{dl), 
where 

fdt. 



G(0 = -^ {P,h-'dwt)-^j^ 

For (0 

Eexp(2G(/)) = Eexp(^-2^ {P,h-^dwt) ~ 



T 

2 



\P\'dt 



T ^ 1/2 



< (Eexp(-4J^ {l3,h-^dwt)-8 \p\^dt) 
x(^Eexp(6 / |/3prft)j 
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C{0,T;H) 





^ d\2 


2 N 1/2 


/ 




dXi) 


JC(0,T;H) 


dXi 



with some constant K^- Hence, 

1 1-^1 ~ -^2 1 1 mr = 2 
1 

< - 

- 2 

< g-7o(fc-m.)T 

for some positive constant 70, since T > is large enough. 

Case 2.S'{{k - 1)T + t, mT) > r + {Ci + l){{k - 1)T + t - mT) for some 
t,i = 1,2. 

This is the direct resuh of Lemma 3.1. In fact, the inequahty 

S\{k - l)T + t, mT) >r+{Ci + l){{k -l)T + t- mT) 
imphes that 

£%{k-l)T+t,mT) > \u{TnT)\'^+Ci{{k-l)T+t-mT) + {r~\u{TnT)\'^) + {k-Tn-l)T. 

Taking 7 = min{l,7o} completes the proof. □ 

Remark 4.2. Lemma 4-i implies that, the probability of the event that the 
orbits in S{m, k — 1), < m < k — 1 enter R{k) U S{k, k) is small. 

To our aim we also need the following lemma. 

Lemma 4.3. Let ^ hold. If {u\u^) G R{k - 1) satisfies \u\{k - 1)T)\ V 
\u\{k-l)T)\ < D, then 

p|(m\m2) g S{k,k)^ >d*>0, 

for some constant d* > depending only on D. 
Proof. From Lemma 3.1, we have 

(17) P|r(t, {k - 1)T) < r + (Ci + l){t -{k- 1)T)} > c* 

for some positive constant c* which depends only on D. For the purpose of 
this lemma, we compare /*(t) with the following standard diffusion process 

(18) l{t) = AI + {0 - u)l + PiF{x,t). 

Here, u is large enough such that (|18|) is linearly dissipative. Let be the 
solutions of (jl8p with the initial value Z*((A; — 1)T). Let Aj(t) and Aj(t) be the 
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distribution of /^.i ^ and restricting on C{{k — l)T,t;Hi), respectively, 

for z = 1,2. Since (fTH| is linearly dissipative, we have 





"ciAi" 


/ 


-dX2- 



for some constant cl > which depends only on the upper bound D. 
By the similar method in the proof of Lemma 4.1, we have 

2 _ 

dXi < C2 i = 1 , 2 



f 


'dXf 




-dXi- 



for some positive constant which depends only on D. Since at this time 
\B''\ = l^r + Pi{r + h'-yi is bounded as we consider the event of (Hlj) (See [TU]). 
By Holder inequality, we can derive 

2 

dXi < cl 



for some positive constant Cg which depends only on D. Then by the Lemma 
C.l of ^7\, we have 

dXi 





\dXii 


/ 


.dX2- 



1 A 



dXr) 



dXo > c*A 



for some constant C4 > depending only on D. Notice that ||Ai — A2I 
1 — ||Ai A A2||^ar5 we have 



Ca. 



I |Al — A2I li.ar < d* — 

Then p|(m^,m^) G S{k^ k)^ >1 — d*. This completes the proof. 
Then, we have 



□ 



Proposition 4.4. Let ^ hold. For u^) G R{k -l)US{k-l,k- I), there 
is a constant < 5 < 1, such that 



p|(u\m2) eR{k)^ <<5p|(m\u2) e R{k-l)u S{k-l,k-l)^ 
Proof. This is the direct result of Lemma 4.3. We omit the proof. 



□ 



By Lemma 4.1 and Proposition 4.4, we have the following important proposi- 
tion, which implies the existence of a unique invariant measure. 

Proposition 4.5. Suppose (OP holds. There is a constant < k, < 1 such that 
for any initial value Mq, Mq G H, 

P| {u\ u^) G R{k) U S{k, A;) I < Kn'' 

where K is a constant depending only on the initial value. 
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Proof. We will consider R{k) and S{k, k) respectively. First, we have 

k-2 

< e-'^'''^^e^'^^¥^{u\u^) eS{m,k-l)^ + 

m=0 

5¥\{v},u^) e S{k-i,k-i)yjR{k-i)^ 

k-2, 



m=0 
m=0 

(e-^^r + 5 + 52)p| (^,1, e - 2) } + 
(252 + e-2^^)p|(ii\ii') eS{k-2,k-2)^ 

where fi, i — 1, 2, 3, increases on k at most polynomially. 
For S{k, k), we have 

F^{u\u^) e S{k,k)^ 

k-2 
m=0 

5P|(ii\ii2) e S{k - l,k - 1)^ + fI^{u\u^) eR{k-l)y 

By the same analysis for R{k), we also have 

This completes the proof. □ 

5. Ergodicity of the Local System 

In this section, we prove the existence and uniqueness of invariant measure, 
and then the ergodicity in the LSSH system. 

Let V{H) be the metric space of all probability measures on {H,B{H)) en- 
dowed with the metric || • ||2 defined by 



sup{^M/i-z/):||/|U<l}, i,,ueV{H) 
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where / is a measurable function on H and 



\l = sup + supj — . . — -.Xy^yeH}. 

X \x-y\ 



It is well known that 1 1 ■ 1 12 generates the weak topology and V{H) is a complete 
space under 1 1 ■ 1 We define the semigroup S{t) on V{H) generated by (jHI)-® 
as 

S{t)^i{T)= [ P{t,x,T)fi{dx) , TeB{H). 
Jh 

A measure /z on is an invariant measure if 



fx{T) = S{t)fi{T) = [ P{t,x,T)fx{dx) , TeB{H). 
Jh 

We restrict S{t) on the set V2{H), which is defined as 

Then we have the following proposition. 

Proposition 5.1. S(t) maps V2{H) into V2{H). 
Proof. For any /i G V2{H) and t > 0, Lemma 3.2 yields 

ljy\'S{t)fi{dy) = l^(^ljy\'P{t,z,dy)y{dz) 



E|-u(t, z)\'^n{dz) 
< e-°* / \z\'^fi{dz) + R 



< oo. 

This completes the proof. □ 
Now we can draw the following main results. 

Theorem 5.2. // (QJj and (0j hold with N large enough, there is a positive 
constant x G (0, 1) such that for any fii, fi2 & V2{H), 

\\S{t)^x, - S{t)fi2\\l < M{fi,,fi2)x' ,t>0. 
Here M {111,112) only depends on fii and fi2- 
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Proof. We follow the proof of [21]. Arbitrarily fixed t > and let k = k{t) be 
the smallest integer such that t < k{t)T, where T is the constant in Sec. 4, 
and let Uq, i = 1,2, be random variables in H with distribution /Xj. Let u^{t) 
and be the weak solutions of (0)-® on [0,kT]. Then we just need to show 
that 

(19) p{t) = F^^\u\t) - u\t)\ > Cie""^*} < 026-"^' 

where Ci, i = 1,2, are positive constants only depending on the initial func- 
tions and (Tj, i = 1, 2, are some positive constants. 

We consider the coupling solution defined in Sec. 4. Define the following event 

[ck] 

Gik) = {iul,ul)e [jSim,k)} 

m=0 

where < c < 1 is some constant. Clearly, 

(20) p{t) <F(G{ky^ +F(^G{k)n{\u^ -u^\ >Cie-"^*}). 
We claim that 

(21) P(^G'(A;)") < Cae-"^* , 

(22) p(g(A;) n {\u^ - > Cie""^*}) = 0. 

First, we prove (24). Since all the orbits in S{m, k) come from R{m)US{m, m), 
we have 

k k 

G{ky C [j S{m,k)UR{k) C [j R{m)US{m,m). 

m=[cfc]+l r?i=[cfc]+l 

Now by Proposition 4.5, we get 

k 

F(G{ky^ < J2 Kte^ <KK^^{k-[ck]) 

m=[cfe]+l 

< KR^ 

for some proper R G (0, 1) and K. Since > t^;, we have 
where cxa = —T^^ In k and C2 = K. 

Next we prove (25). It is enough to prove that if {v},v?) G G{k) and Ci is 
large enough 

(23) - < Cie-"' 



(72 < 
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Indeed, by the definition of G{k) for {u^,u^) G G{k), there is an integer /, 
< / < [cfc], such that u^) G S{1, k). Therefore, the relations (|I2I)-(|I1D are 
satisfied. Then as Proposition 2.1, for u = Ui — U2, we have 



\u 



{t)\ = \pit)\<2de~'^'~'^\ 



Notice that IT < ckT < c{t + T), then t - IT > {1 - c)t - cT. Hence, 
\u{t)\ < 2rfe^^^e-^(^-'=)*. Thus (25) holds with Ci = 2rfe^^^ and ai = A(l - c). 
Finally, taking M{fii,fi2) = C2 and x = C^'^ , we obtain the conclusion and 
the proof is completed. □ 

Corollary 5.3. // enough modes are forced and the condition (OJ) holds, then 
the LSSH system has a unique invariant measure /ig, such that for any u & H 
and t > 0, 

\\P{t,u,■)-^l^\\l<M{\u\^)x\ 
where M(|up) is a positive constant and x is as in Theorem 5.2. 

6. Ergodicity of the Nonlocal System 
In this section, we turn to the following NLSSH system on the real line 

(24) ut- Qu + {1 + d^^fu + uG*u^ = F{x, t) , 

(25) m(0) = Mo, 

(26) u{t,x) = u{t, x + 2n), 

where 

G*u^= / G{x-Ou\^,t)d^. 
Jo 

We always assume that 

(27) bi^O , l<i<N<N" 



N" 

2 



(28) 5o = ^6; 



i=l 



First, we consider the following case: for every x G M 
(29) < 6 < G{x) < a, 

where a and b are the positive constants (i.e., G > is a positive kernel). Then 
we have the same energy estimation as LSSH system. 

Lemma 6.1. There is a positive constant Ci, such that for any r > 0, 
p||M(t)P + ^ \u{s)\^ds< |m(0)|2 + Cit + r ,/ort > 0} > 1 - e'^ 
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Proof. We only estimate the nonlinear term. In fact, 

{uG*u\u) = {G*u\u^) > b\u\\ 

Set Ci = max^, |2(^ + 6^^^ + - + 5o|, where 
(30) 6max = max6j. 

i 

By using the same analysis as in Lemma 3.1, we can finish the proof. □ 

Note that, for the nonlocal nonlinearity, only the part of results in Proposition 
2.1 still holds. But this does not affect the coupling result. Taking the same 
decomposition of NLSSH system as one of LSSH system, we have 

Proposition 6.2. For any u^, -Ug G H, if ll^^i^ = ^ and , v? satisfy |T^, 
then there exists a constant A > such that 

W-u'^\ < e-'^'^'-'^'^^\u\mT) - u'^{mT)\ , t e [mT,kT]. 

Proof. Since ll^k — ^mk^ '^^^J consider the high mode. Let p = — . 
Then, p satisfies 

p = Ap + Qp- Phiu^G * {u^f - u^g * (u^f). 
Taking the scalar product of the above equation with p in Hh and notice that 
{u^G*{u^y -u^G*{uy,p) 
= {u^G * {u^f - u^G * {u^f + u^G * {u^f - v?G * {v?f, p) 
= {u'G * {{u'r - {u^ f),p) + {{u^- u')G * {uy, p) 
> {u^G*iu^ -u^){u' + u^),p), 

we have 

dt\p\' < {af, + g)\p\' + \{u'G*{u'-u'){u'+u'),p)\, 
< {af, + g)\p\^ + a{^\u'\' + \u'\')\p\^ 
where is the Nth eigenvalue of A. Then 
br < \p{mT)\'^exp\aj^ + g + a [ Au^l"^ + \u^\'^)ds] 

< |p(mr) p exp {[a^ + g+ /^^l^'l' + l^'l')^^) " "^^) } 

< \p{mT) 1^ exp I (a^ + g + ^aCi + a) (t - mT) | , 

since we can always take r < T. If we choose so large that a^+^+|aC*i+a < 
0, then 

|^|2 < e-^(*-'"^)|p(mT)| , t e [mT,kT], 
for some constant A. This completes the proof. □ 
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Remark 6.3. The above proposition is weaker than Proposition 2.1, hut it is 

enough for the coupling result in Lemma 4-1- 

We have to give more estimation for the couphng result. For Lemma 4.1, we 
have to estimate \B{1, h},, hl)\ = \Piu^G * {u^ - Piu^G * {u'^Y\. In fact, 

\u'G * iu'y - u^G * (u^l 
= \{l + h')G * {I + h')^ - (/ + h^)G * {I + h^)^\ 

< \lG*[{l + h'y -{l + hy]\ 

+ \h^G * {f + {h^f + 2lh^) - h^G * {f + {h'^f + 2//i2)| 

< \21G * \l{h^ - + \IG * \{h^ - h^){h^ + 

+ \{h^ - h^)G * P\ + \h^G * {hy - h^G * {h^l 
+2\h^G * {Ih^) - h^G * {lh^)\. 

From ()29|) we have 

(31) \21G * [l{h^ - h^)]\ <2a\h' -h^\\l\''. 



(32) \IG * [{h^ - h^){h' + < a\h' - h^\\h' + h^^ 



(33) \{h' -h^)G*l^\<a\h' -h^\\l^\. 

(34) \h'G * {h'f - h^G * {h^l 

(35) < \h'G*{h'f -h^G*{h^f\ + \h^G*{h'f ~h^G*{hy\ 

(36) < a\h^ - h^\\h^\^ + a\h^ - h^\\h' + h^\\h% 

Similarly, we have 

(37) 2\h'G * (Ih^) - h^G * {lh^)\ < 2a\h^ - h''\\h^ + h^\\l\. 

Then, from (33)-(39) and |mP increases polynomially, we can derive the cou- 
pling result in Lemma 4.1 by Proposition 6.2. 

For Lemma 4.3, we have to estimate \B\ = \vl + PiuG * (Z + hY\. Since 

(38) \uG * u^\ < a\u\^ , 
\B\is bounded. 

With all the above analysis, we draw the following conclusion. 



Theorem 6.4. // |F?D and / f^) hold with N large enough, the NLSSH sys- 
tem with the positive kernel h2f^) has a unique invariant measure, and it is 
exponentially attracted in V2{H). 
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The above analysis for the positive kernel ()29|1 does not work for the non- 
negative kernels. But we can still work for a special nonnegative kernel. See 



Define 

(39) J(x) 
where 



cexp(-Yz^) if X < 1, 
ifx>l. 



c={ I exp{-^-—^)dx 



We further define that for 5 > 0, 



Jsix) = 6-'Jij). 



Let Co(/) be the space of continuous functions with compact support in I. It 
is known that 



(40) 11^-5* /-/lbo(/)^0 as 5^0, 

for any / G Co(/). Thus for any e > 0, there is a 5o = (5o(e) > 0, such that 

(41) J5o*f>f-e. 



We consider a special kernel G{x) = Js^ix), with Sq = 6o{e) as in (jH]), for the 
NLSSH ©. Then, G satisfies 

(42) < G < |. 

Note that the lower bound of G is 6 = 0, so some estimation above does not 
apply. But due to PT|) we can still get the similar energy estimates. And 
notice that we only use the lower bound in the energy estimates, so we can 
have the same ergodic result. 

From [IH], or we know that the problem has a unique 

solution u{t,x,uo) e C{[0,T]] H) n L'^{{0,T), H^{I)) and 

u{t,x,Uo) G Go(/),Vt > to > 0. 

Thus 1)411) implies that 

(43) G{x) * > - e. 
Then the estimates of Lemma 6.1 still holds with 

(44) Ci = max |2(/i + + ^ + e)x^ - ^x" + 5o}. 

The other estimates does not depends on the lower bound of G. Then we can 
derive the following result. 
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Theorem 6.5. // jFTp and / T^) hold with N large enough, the NLSSH system 
with the nonnegative kernel Js^ix) has a unique invariant measure, and it is 
exponentially attracted in V2{H). 

7. Dynamical Difference between Local and Nonlocal Systems 

With the resuhs on ergodicity of the local and nonlocal stochastic Swift- 
Hohenberg equations in previous sections, we can now compare some dynam- 
ical behavior. 

Prom Proposition 2.1, we see that the local stochastic Swift-Hohenberg system 
is ergodic provided aN + Q < 0. Note that aN = —{N"^ — 1)^. So we only need 
the number N of randomly forced modes to satisfy 



Note that the number N does not depend on the random forcing term, i.e., it 
does not depend on the coefficients b'^s in the random forcing in (jSj). 

However, from Proposition 6.2 the nonlocal stochastic Swift-Hohenberg system 
with positive kernel (0 < 6 < G{x) < a) is ergodic provided that satisfies 

aj^ + g+ laCi + a < 0, where ^1 = 50+ (^±^:^£±il!, with Bq defined in ^ 
and hmax defined in ()30|) . Thus we only need the number of randomly forced 
modes to satisfy 



Clearly, we have the following comparison between the local and nonlocal sto- 
chastic Swift-Hohenberg models: (i) The number of Pourier modes to be ran- 
domly excited for ensuring ergodicity of the local stochastic Swift-Hohenberg 
system depends only on the Rayleigh number through > 0, which measures 
the difference of the Rayleigh number from its critical convection onset value; 
(ii) The number of Pourier modes to be randomly excited for ensuring ergod- 
icity of the nonlocal stochastic Swift-Hohenberg system depends on the bound 
of the kernel G in the nonlocal term, and the random term itself, as well as on 
the Rayleigh number. 





Por the nonnegative kernel G{x) with upper bound a, we only need 
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